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ABSTRACT 

We present a formalism of the dynamics of internal shocks in relativistic jets where 
the source has a time-dependent injection velocity and mass-loss rate. The variation 
of the injection velocity produces a two-shock wave structure, the working surface, 
that moves along the jet. This new formalism takes into account the fact that momen- 
tum conservation is not valid for relativistic flows where the relativistic mass lost by 
radiation must be taken into account, in contrast to the classic regime. We find ana- 
lytic solutions for the working surface velocity and radiated energy for the particular 
case of a step function variability of the injection parameters. We model two cases: a 
pulse of fast material and a pulse of slow material (with respect to the mean flow). 
Applying these models to gamma ray burst light curves, one can determine the ratio 
of the Lorentz factors 72/71 and the ratio of the mass-loss rates m^jrax of the up- 
stream and downstream flows. As an example, we apply this model to the sources GRB 
080413B and GRB 070318 and find the values of these ratios. Assuming a Lorentz fac- 
tor 71 = 100, we further estimate jet mass-loss rates between m,\ ~ 10~ 5 — 1 M yr~L 
We also calculate the fraction of the injected mass lost by radiation. For GRB 070318 
this fraction is ~ 7 %. In contrast, for GRB 080413B this fraction is larger than 50%; 
in this case radiation losses clearly affect the dynamics of the internal shocks. 
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1 INTRODUCTION 

Collimated outflows (with jet-like geometry) moving at rel- 
ativistic speeds are characteristic of active galactic nuclei. It 
is commonly accepted that an extragalactic jet is produced 
in the neighborhood of a massive black hole in the center of 
an active galaxy (e.g. Rees 1984; Istomin 2010). These rela- 
tivistic jets are subject to the development of shock waves. 
Rees (1978) proposed that the observed knots in the ex- 
tragalactic jet M87 correspond to the locations of internal 
shocks which arise owing to variations in the outflow velocity 
of a beam generated in the nucleus. Later, Rees & Meszaros 
(1994) pointed out that fluctuations of the Lorentz factor 
around its mean value in a relativistic outflow, that give rise 
to internal shocks, can dissipate a substantial fraction of the 
outflow energy into non-thermal radiation. They proposed 
that this mechanism is operating in the so-called gamma-ray 
bursts (GRBs). 
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Several authors have studied internal shocks in ultra- 
relativistic outflows to explain the observed variability 
of GRBs (e.g., Mochkovitch, Maitia & Marques 1995; 
Kobayashi, Piran & Sari 1997; Daigne & Mochkovitch 1998; 
2000). In these models the flow is represented by a succes- 
sion of shells with different values of the Lorentz factor. This 
models reproduce the burst profiles and their short-time 
scale variability. Kobayashi et al. pointed out that variations 
of the relativistic flow velocity are strongly correlated with 
the temporal variations observed in the GRBs. Daigne & 
Mochkovitch (1998) studied the detailed radiation processes 
to calculate the fraction of the kinetic energy dissipated in 
the shocks that can be emitted in the form of gamma rays 
and obtained that the total efficiency is of the order of only 
a few percent. In addition, Spada et al. (2001) proposed that 
the internal shock scenario can also be used for blazars. For 
comprehensive reviews of the physical processes and obser- 
vations of GRBs see, e.g., Meszaros (2002); Piran (2004) ; 
and Gehrels et al. (2009). 

Using mass and momentum conservation, Canto et al. 
(2000) solved the dynamics of internal shocks of non rela- 
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tivistic jets with time dependent injection velocity and mass- 
loss rate. Mendoza et al. (2009) used this momentum con- 
serving formalism in the case of relativistic jets and com- 
pared with observed light curves of GRBs assuming a sinu- 
soidal velocity variation. However, momentum is not con- 
served in relativistic flows because radiative losses change 
the relativistic mass since the Lorentz factor decreases when 
energy is radiated away. In this paper we present a new for- 
malism that describes the dynamics of internal shocks in a 
relativistic jet taking into account the momentum change by 
radiation. In particular, we study the dynamics of internal 
shocks for the case where the injection velocity and mass- 
loss rate are both step functions of time, when a fast wind 
reaches a slower flow. For this type of variability, we find an- 
alytic solutions for the dynamical evolution and luminosity 
of the shocks, which are implicit functions of time. 

The organization of the paper is as follows. In §2 we 
discuss the relevant relativistic equations. In §3 we present 
the model, and investigate the dynamical evolution of the 
internal WS. The luminosities predicted by our model from 
relativistic jets are given in §4. A comparison between our 
analytic solutions and observations of extragalactic gamma 
ray bursts are presented in §5. In §6 we evaluate the mass 
lost by radiation in the WS. Finally, in §7 we summarize our 
conclusions. 



2 RELATIVISTIC EQUATIONS FOR 
INTERNAL SHOCKS 

For a free-streaming flow, the non-dimensional velocity at a 
distance x from the source, at time t measured at the source 
reference frame, is given by 

x 



WS 



P(x,t) =/3(0,r) = 



c(t - 



(1) 



where the injection point is at x — 0, and r is the time at 
which the flow was ejected, and, as usual, /? = v/c, where 
v and c are the speeds of of the flow and of light, respec- 
tively . If the flow velocity at the injection point increases, 
the fast material will reach the slower material and form a 
working surface (WS) bounded by two shock fronts (Raga 
et al. 1990). These WS structures are the so-called "internal 
shocks" . Figure[T]shows a schematic diagram of a WS formed 
when a fast upstream material with velocity fa reaches a 
slower downstream material with velocity fa . The WS moves 
with velocity {3 WS intermediate between fa and fa. 

The internal WS forms at the distance Xi from the 
source given by 



d/3{r)/dT 



(2) 



where we wrote for simplicity, fi{r) = /3(0,r), and the min- 
imum is taken over the time interval where the velocity 
increases. The WS is formed at time U from the material 
ejected at time n, both values given by equation Then, 
one can obtain the time, tt, at which the WS is formed from 
equation fTJ). 

Consider a relativistic jet with time dependent injec- 
tion velocity, /3(r), and mass-loss rate, rh(r). As discussed 

1 First, tj is obtained by minimizing the RHS of this equation. 
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Figure 1. Schematic diagram showing a working surface formed 
by the interaction between two relativistic flows. The upstream 
and downstream flow velocities are fa an d fa (with fa > fa), 
respectively. The working surface moves with an intermediate ve- 
locity p wa . 



above, a WS is formed when fast material overtakes slow 
material. This WS travels downstream the jet flow with a 
velocity /3 ws (t), where the time t is measured in the source 
reference frame. The slow and the fast material just enter- 
ing the WS at time t were ejected at times n and T2, re- 
spectively, with corresponding downstream and upstream 
velocities fa — /3(ri) and fa = P{t2), and mass-loss rates, 
rhi = rin(ji) and rii2 = Jtc(ts). The time dependence of the 
velocities and mass-loss rates is given through the time de- 
pendences Ti(t), and T%(t). 

In Appendix A we discuss a simple example of the in- 
elastic collision of 3 relativistic particles that radiate energy 
after they collide. This example shows that the momentum 
of the final particle is not conserved because its relativistic 
mass changes when energy is radiated away. Therefore, we 
introduce below the energy and momentum equations that 
take into account the energy lost by radiation. 

The total energy dissipated by the flow interaction E(t) 
is given by the difference between the total energy injected 
into the WS and the energy carried by the WS at the instant 



E(t) 



m(r)7(r)dr - m ws (t)7 ws (t), 



(3) 



where 7ws(i) = Vv^ ~~ Pws(t) 2 , and the rest mass injected 
into the WS is 



m wa (£) = / m(r)dr. 



(4) 



We assume that dissipated energy E(t) is completely 
radiated away (see eq. |A9] 1: none of this energy is stored in 
internal degrees of freedom. Therefore, the luminosity (the 
radiated energy per unit time in the WS) is given by the 
time derivative L r (t) — dE(t)/dt. Then, the dynamics of 
the WS is described by the energy equation, 

j t (m ws (t.) 7ws (t)) = j t f m(r) 7 (r)dr - (5) 



2 For simplicity, the internal energy of the particles that enter 
the WS is ignored. 
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Figure 2. Injection velocity /3 as function of time r: I) The initial 
velocity /?i suddenly increases to /3 2 at a time r = 0, for a finite 
time interval, At, and then instantly returns back to its original 
value. II) The initial velocity /9 2 instantly decreases to /?i at a 
time r = —At, and at t = 0, the faster flow starts to be injected 
again. 

obtained from the derivative of equation (|3}, and the mo- 
mentum equation is given by 

m(r)7(r)/3(T)dr 

' ti 
L r (t) 



(m ws (t)7„ s (f)/3 ws (t)) = 



-Mt), 



(6) 



where the last term in the RHS is the momentum change due 
to the relativistic mass lost by radiation. Combining equa- 
tions ([5]) and ([6]) one obtains the equation for the velocity 
of the WS, 



3 (t)7ws(£) 



dt 



d f 

= -T. \ m( T h( T )l 3 ( T ) dT 

J T\ 



-/M*)| 



-m(n)7(n)[*)-fc(i)]|, 



= m(T 2 )7(T 2 ) [/3(r 2 ) 



that does not depend explicitly on the radiated energy, L r . 

Given a time variability of the functions /3(r) and m(r), 
one can find the WS velocity, /3 ws (i), from equation © to- 
gether with equation (fTJ) , the latter giving the relation be- 
tween ti,2 and t. We will show below how these equations 
can be solved analytically for the case of step function vari- 
ability of the injection parameters. 



3 STEP FUNCTION VARIABILITY OF THE 
INJECTION PARAMETERS 

Consider a step function variability for the injection velocity 
and mass-loss rate such that, for r < 0, a slow flow has 
ft = p\, and m = rhi; and for r ^ 0, a fast flow has /3 = /3 2 
and m = m 2 . Figure[2]shows two cases: I) when fast material 
is injected at r > for a finite time, At; and II) when slow 
material is injected at r < 0, for a finite time interval, At. 
In both cases, when the fast flow reaches the slow flow, the 
WS is formed instantaneously at the injection point: x — 0, 
t = 0, and t = 0. 

As we will show below, the dynamical evolution of the 
WS goes through 2 stages. In the first stage, the WS is fed by 
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Figure 3. Qualitative behavior of the WS velocity as a function 
of normalized time t/ti' 11 , for Case I and Case II, respectively. 
See text for description of this figure. 



both the slow and fast flows and it moves at a constant ve- 
locity, intermediate between the fast and slow speeds. The 
second stage begins when one of the flows has been com- 
pletely incorporated into the WS; then, the WS accelerates 
or decelerates depending on which flow (fast or slow) contin- 
ues to feed the remaining shock. Asymptotically in time, the 
speed of the WS in the second stage tends to the velocity of 
the remaining flow. 

Figure [3] shows the qualitative behavior of the WS ve- 
locity as a function of normalized time for Case I (pulse of 
fast material) and Case II (pulse of slow material). In Case 
I, the time is normalized to the critical time t], (eq. P3]). ^ n 
Case II, the time is normalized to the critical time t 1 / (eq. 
[26]). In both cases, the constant velocity phase ends when 
t/#" = l. 

Consider the material at time t that enters the WS, lo- 
cated at the position x wa (t), through both shocks. The slow 
and fast material were ejected at times n and r 2 , respec- 
tively, such that, according to equation JT} 



Tl,2 



t 



Pi 



and 



dt 



= 1 



(8) 



where fi,2 = ti, 2 /At, t = t/Ar, and £ ws = i ws /(cAt). 

Now, for a step function variability, the rest mass of the 
WS, given by equation @, is 



m wa (i) = 



l-T 2 

rhi(T)dT + / rh 2 (r)dr 
n -'o 
■rom(*) + ra 2 T 2 (i). 



(9) 



Note that n < and t 2 > 0, and for now on, for sim- 
plicity, we will drop the t dependence of the functions. 
The velocity of the WS is obtained from equation 

(/3 2 - /? ws ) 2 



dt 



m 2 7 2 - 



m i7i 7, , (10) 

Pi 



where we used equation j5J. This equation has the constant 
velocity, /3 ws o, solution 0, such that d/3 wa o/dt = 0, then, 

A(/3 2 -/3 W80 )=±(/3 w8 o-/3i), (11) 

with 



A = 



rh 2 7 2/ 9i 
mi7ip2 



(12) 



3 This solution exists because at t = 0, the initial condition /3 W sO 
makes the RHS of the equation equal zero, since m ws (0) = 0. 
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where we defined the velocity ratio a — P2/P1, the mass- 
loss rate ratio b = mil mi, and gamma ratio r — 72/71. 
The correct solution corresponds to the + sign Q, with the 
ordering Pi < /? ws < P2, where the constant WS velocity is 
given by 

A/? 2 + Pi 



fiwsO — 



7u)s0 



A + l 



v /(A/ 7 2) 2 + 2A(l-/3 1 /3 2 ) + l/7i 2 



(13) 



This velocity corresponds to the first stage of the evolution 
of the WS, when it is fed (and bounded) by two shocks. 

Now we will discuss the evolution of the WS in the sec- 
ond stage of the two different cases I and II shown in Fig- 
ure [2] Even though the formalism is the same, the resulting 
equations are different for each case. Thus, for clarity, we 
separate them in the next two subsections. 

3.1 Case I: Decelerating WS 

The constant velocity phase ends when T2 = At, i.e., when 
the fast material is completely incorporated into the WS. 
This happens at a critical time obtained by substituting the 
position of the WS, x wa o = p ws ot, into equation (JSJl, 

_ a (A + 1) 

If. — 



a- 1 ' 

which corresponds to ejection time 



TV. 



Pi 



-aA. 



(14) 



(15) 



In this second stage, t>V c , the first term of equation 
is because dri/dt = 0. Furthermore, substituting T2 = At 
in equation @, the rest mass of the WS is 



= tu^At — m\T\ = rhi At (b — fx) 



(16) 



Collecting these results we write the WS velocity in equation 
as function of fi as 



71 A) 



(17) 



dfi 7 {b - h) 

For /3 WS 7^ Pi , this equation^ can be integrated by separation 
of variables as 

2 (71 - PwsPlJl + 1/7 

(p—pT) = c ( b --0, (is) 

that is a quadratic equation for /3 W s(ti), and has the solution 

, PiT 2 +4T + 4Pi 

Pw S (n) = „ . - ^ . and 



T 2 +APxT + A 



7w a (n) = 71 



T 2 +4/3iT+4 



7-2-4 

where we defined the time function 
C'ib-fx) 



T(n) 



71 



(19) 



(20) 



4 The — sign in equation lilt is unphysical because, in this solu- 
tion, /3 WB o < 01 , that implies that there is no downstream shock. 

5 Equation 1171 has a constant velocity solution, /3 W s = Pi, that is 
trivial because there is no shock. 



that is an increasing function of |fi|. From equation (|18|) one 
can show that T(fi) > 2. The constant C 1 is obtained by 
matching the solution /3 ws (fi jC ) = P wa o, where fi, c is given 
by equation (I15|) ; thus, from equation (|18[) one gets 



(21) 



r i _ 2 (71 - PwsoPi-yi + Vtwso) 
(b + aX)(p ws o-Pi) ' 

To find a relation between fi and t, we take the derivative 
of the time function T (eq. |20]). 

dT C 1 dfi _ 4C 7 T 
~dl ~ ~~~dl ~ Pi-il (7^+4/?iT + 4)' 



(22) 



where we used equations © and ()19p . Again, by separation 
of variables, this equation can be integrated as 

2 Pill 



(23) 



where the constant D 1 is obtained evaluating this expression 
at tr c and fi, c , and is given by, 

2 A7i ( a - 1) 

where the critical time function is 

T c = T(fi, c ) = C' (b + aA) 



(24) 



7i 



(25) 



Therefore, given fi one can evaluate the time function 
T to obtain /3 W8 (fi) from equations (fT9)l and (f20)l . and one 
can obtain the WS velocity as a function of time t, P ws (t), 
in a tabular form, using equation (|23|) . Finally, the position 
of the WS, Xvrs(t), is given by equation ((8|. 



3.2 Case II: Accelerating WS 

In this case, the constant velocity phase ends when n = 
— At, i.e., when the slow material is completely incorporated 
into the WS. From equation (JSJ, the critical time is 



~ t u = (A + l 



A (a — 1) ' 
corresponding to the ejection time 

PvisO \ i 



T2,c = t'J I 1 



Pi 



aX 



(26) 



(27) 



For t > t c , one has dri/dt = 0. Also, from equation 
((9|, the rest mass is 



rii2T2 + mi At — rili At (&T2 + 1) . 



(28) 



Then, we write the WS velocity in equation as function 
of T2 as 



dP^ 



72 K A -M 



(29) 



df 2 7ws (bf 2 + 1) 

For p wB 7^ P2, this equation can be integrated by sep- 
aration of variables as 



2( 72 - /?wsA;72 + l/ 7 ws) n u th ~ , 

= C (&T2 + 1) - 



(30) 



B Equation 1291 has a constant velocity solution, /3 WS = P2 , which 
is trivial because there is no shock. 
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which is a quadratic equation for /3 W s(t 1 2), and has the solu- 
tion 

ftT 2 - 4T + 4/3 2 



/3ws(t2) = 



r 2 -4/3 2 r + 4 



and 



7ws(T2) = 72 



r 2 -4/? 2 r+4 



7-2-4 

where the time function is 

C"(bT2 + 1) 



72 



(31) 



(32) 



that is an increasing function of T2. As in the previous case, 
one can show that 7"(t2) > 2. The constant C is obtained 
by matching the solution /3 ws (t2, c ) = ftso, where fa,c is 
given by equation (f2T)l . Substituting these results into equa- 
tion (|30|l one gets 

_ 2a A (72 - /3 W5 o/3 2 72 + 1/7wbo) , , 

(6 + aA)(ft.-/3 ws0 ) ' 1 J 

We take the derivative of the time function T (eq. [32]), 



dT C n b df 2 AC" 



dt 



T 



72 dt /?272 (T 2 -4/? 2 T + 4) : 



(34) 



where we used equations (JSJ) and (|31|l . By separation of vari- 
ables, this equation can be integrated as 



7-2 4pii 

-7; 4/?2T + 41nT= -3 — jt + D n , 

2 ft 72 



(35) 



where the constant D /7 is obtained evaluating this expres- 
sion at and T2, c , and is given by, 



2 /?275 A (a - 1) 

where the critical time function is 
// (& + «A) 



7^ = T(f 2 , c 



aA72 



(36) 



(37) 



Using equations (|31|) . (|32[) . and (|35[) one can proceed as 
in Case I to obtain /3 ws (t) and x ws (t) as functions of t in a 
tabular form. 



4 LUMINOSITIES 

Using a step function variability of the injection parameters 
in equation ((5j) , the luminosity L r of the WS is given by 



L r (t) 



dfo 



dri 



c 2 dt dt 

, / . df 1 . df2 \ m ws d7ws 

+ mi — - =- — m 2 — - 7„ s 

V dt / Ar dt 

. - s df 2 . . , x dn 

= m 2 (72 - 7ws) +mi (7w S - 71) -1 f 
dt dt 



m„ s d7 W c 
Ar dt 



(38) 



where we have used equation ([9]) for the rest mass m ws - 

Only a fraction of the energy radiated in internal shocks 
will be emitted as gamma rays; this fraction is low, e ~ 0.01 
(e.g., Daigne & Mochkovitch 1998). Here we will assume that 
a constant fraction e of the luminosity will go into gamma 
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Figure 4. Normalized luminosity L/L°, as function of the nor- 
malized time t/t c for Case I and Case II. For these models, we 
assumed 71 = 100 and 72 = 200. We also assumed A = 1, thus, 
the critical times are the same for both cases. 



ray radiation, Lqrb = eL r (t), and express this GRB lumi- 
nosity in non dimensional form as 



L 



Lgrb 
ec 2 7imi 

\2 (72 - 7ws) df 2 (7„ s - 71) dfi 
— &A — -)- 



72 dt 
m ws d7ws 
mi7iAr dt 



7i 



dt 



(39) 



In the constant velocity phase, i < t c ' , the luminosity is 
given by 



^0 ,2 (72 - Jwao) ! ftsO 



72 \ ft 

(-fwsO — 7l) / ftsO 



71 



ftsO 



- 1 



ft 



(72 - 7»s0) {"fwsO - 7l) 



72 



71 



(40) 



where we substituted dfi/dt = 1 — /3 ws //3i and dr^ldt = 
1 — /3 ws /ft from equation (JSJl. 

As mentioned in §3.1, in Case I, once all the fast mate- 
rial has been been completely incorporated into the working 
surface, it is decelerated by the slow downstream flow. This 
second stage starts at a time ti given by equation (|14[l . For 
this case, on has dr^fdt — in equation (|38p ; thus, one can 
write the luminosity as, 



£J _ (7W8 - 7l) ( ft»«0 _ 1 \ _ ^ d7w B 

rhi7iAr dt 

2 (7ws - 71) 



7i v ft 

(ft™ -ft) 



ft 



ft s (ft™ -ft) 7 



7i 



(41) 



where we used that d7 ws /dt = /3 ws 7^ s d/3 ws /dt, and, from 
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equation (0, 

df3 ws _ mi7iAr (/3i - (5 WS ) 2 



dt 



^ws7ws 



Pi 



(42) 



On the other hand, in Case II, once the slow material 
has been completely incorporated at FJ given by equation 
(|26p . dfi/dt = 0, and the luminosity is given by 



fll ,2 (72 — Jwso) I PwaO 

72 

\ 2 a ^ 2 - ^ ws ^ x 



mi7iAr dt 



fa 

"/3ws (fa - fas) 7 + 

where, from equation (J7J), 

d^ wa _ 771272 Ar (#2 - /9ws) 2 



(72 - 7ws) 



72 



(43) 



(44) 



dt m 7ws fa 

In the ultra-relativistic (UR) limit, 7 > 1, the expres- 
sions for the luminosities are simplified as shown in Ap- 
pendix B. In the following section we will apply these equa- 
tions for the GRB luminosities to describe the light curves 
of two observed sources. 



5 PREDICTED FLUXES 

Figure [4] shows the model luminosity for Case I and Case II, 
normalized to the luminosity in the constant velocity phase, 
L 1 ' 11 /L° , as a function of normalized time, t/tjj.' 1 . The mod- 
els have 71 = 100, 72 = 200, that correspond to the UR 
limit. Also, we choose A ~ (br) 1 ^ 2 — 1, that implies a con- 
stant energy injection rate rhiji = 771272. In this case, the 
critical times (eqs. [14] and [26]) are equal. For this reason, 
we drop the superscripts in the following discussion. As dis- 
cussed above, for t < t c , the WS is bounded by 2 shocks and 
moves at a constant speed. Then, the relative velocity be- 
tween the incorporated material and the WS and, therefore, 
the luminosity are constant. For t ^ t c , in both cases one 
shock disappears and the relative velocities between the WS 
and the new material that is incorporated decreases with 
time; therefore, the luminosity, diminishes with time. 

Appendix \C\ and [D] show that, in the UR approxima- 
tion, one can fit a power-law to the wing of the GRB light 
curve, and obtain the critical time function T c , the gamma 
ratio, r — 72/71, and the mass-loss rate ratio, b = 7712/7711. 
These quantities can be obtained without any further as- 
sumptions. As an example, Figure [5] shows a model fit to the 
observed sources GRB 08413B and GRB 070318. The fluxes 
were taken from Mendoza et al. (2009) that correspond to 
observations between 15 and 150 keV with the Burst Alert 
Telescope on board the SWIFT satellite. We first fit the ob- 
served flux density F directly and later discuss its relation 
to the distance and luminosity of the GRB. We follow the 
procedure described in Appendix [D] and choose the value 
of the luminosity (or flux density) in the constant velocity 
phase, the time to of the beginning of the velocity pulse, and 
time t e = t c + to of the end of the constant velocity phase, 
where t c is the critical time of the model. We also choose the 
flux at the constant velocity phase, _F°. The wing light curve 
is then fit by a power- law F = B(t — to) a . As discussed in 



the Appendix C, the value of the slope a, determines which 
case (I or II) applies. 

Table 1 shows the model parameters that fit the light 
curves of both GRBs. The name of the GRB is indicated in 
column 1; the applied model (Case I or Case II) is shown 
in column 2; column 3 shows the flux in the constant ve- 
locity phase; column 4 and column 5 show the time of the 
beginning of the pulse, to, and the time of the end of the 
constant velocity phase, t e ; column 6 and column 7 show 
the coefficient, B and the exponent a of the power-law fit 
of the GRB light curve wing; column 8 gives the inferred 
value of the critical time function, T c \ column 9 gives the 
jump of the flux (or luminosity) at the end of the constant 
velocity phase, J(t e ) defined by equations (IB9|) and (|B10|) ; 
finally, column 10 and column 11 give the inferred values of 
the gamma ratio, r and the mass-loss rate ratio, 6, obtained 
directly form the fit to the wind of the light curve. 

In our model, the WS approaches the observer at rel- 
ativistic speeds. As discussed in Appendix [E] the observed 
bolometric flux density of an approaching relativistic jet is 
increased with respect to an emitter at rest by a factor 8^ s , 
where 5 WS = l/[~y vs (l — p vs cos 8)] and 9 is the angle between 
the direction of the relativistic jet and the observer. Lind & 
Blandford (1985) obtained an amplification by a factor of 
8^ a between the observed flux at frequency v and the emit- 
ted luminosity at frequency v 1 = v/8 ws . An extra factor of 
5 WS is obtained when one integrates the frequency to get the 
observed bolometric flux (or the observed flux in a frequency 
range) in terms of the emitted bolometric luminosity. Thus, 
in our model the observed gamma ray flux in the constant 
velocity phase is given by Fq RB = 8% s oL^j R ec?^iirni/4,-KD 2 , 
where L% R is the normalized luminosity in equation (|B5|) . 

In order to solve for the mass-loss rate, we assume that 
the relativistic jet is seen almost along the jet axis, i.e., 
cos9 ~ 1, thus, <5 W sO ~ 27 W so. Also, assuming 71 = 100, both 
L\, R and the Lorentz factor 7 ws o can be obtained from the 
model parameters in Table 1. In particular, we obtain large 
Lorentz factors for the WS, 7 ws o ~ 802 for GRB 080413B, 
and 7 ws0 ~ 120 for GRB 070318 from equation fB2]). Fur- 
thermore, we estimate the distance to the GRBs from their 
redshift. The source GRB 070318 has an estimated redshift 
z = 0.836 ( Jaunsen et al. 2007) and GRB 080413B has a red- 
shift z = 1.1 (Vreeswijk et al. 2008). Assuming a dark energy 
density Qa = 0.7, a matter density f2 m = 0.3, and a Hub- 
ble constant H — 70kms _1 Mpc _1 , the luminosity distance 
of GRB 070318 is D = 5300 Mpc and of GRB 080413B is 
D — 7440 Mpc. With all these ingredients, we can now solve 
for the mass-loss rate, and obtain 7771 ~ 9.7 x 10 -6 Mq yr _1 
for GRB 080413B, and mi ~ l^Moyr -1 for GRB 070318. 
Because the emitting WS moves with relativistic speed, the 
jet mass-loss rates mi required to produce the gamma ray 
flux observed at the Earth are much smaller than those 
obtained, for example, in nucleosynthesis models of wind- 
driven supernovae and collapsar models (e.g., Fujimoto et 
al. 2008; Maeda & Toming 2009). 

Note that in the model the time is measured in the 
frame of reference of the jet source, where the evolution 
timescales for the WS are of the order of ~ 10 5 — 10 6 s. 
Instead, the time of the observations is measured at the ob- 
server's frame of reference and is only of the order of tens 
of seconds if the relativistic jet is seen almost along the 
jet axis. This happens because the arrival time is At a = 




GRB 070318 
CASE II 




240 250 260 270 280 

t(s) 



Figure 5. Left panel: GRB 080413B source with a decelerating WS model (thick solid line). Right panel: GRB 070318 source with an 
accelerating WS model (thick solid line). The thin lines in each panel show the duration of the constant velocity phase. The observational 
data of both sources was taken from Mendoza et al. (2009). 



Table 1. Model Parameters 



GRB 


Case 


p0 a 
GRB — 

erg cm _2 s _1 


to 

(s) 


te 
00 


B±AB 
(10 _7 erg cm" 2 s" 1_ °) 


a ± Aa 




J 


r 

(72) 


b 

(») 

v mi ' 


080413B 


I 


1.08 X 10" 6 


238.50 


240.01 


15.49 ±3.00 


-2.03 ±0.30 


2.57 


0.38 


22.49 


233.73 


070318 


11 


1.25 x 10~ 7 


239.00 


245.90 


8.08 ±0.50 


-1.08 ±0.05 


5.35 


0.20 


2.64 


0.12 



Flux density in the constant velocity phase, F£U 



ec 2 71 mi / 4-7T D 2 



At em /(7 wa J ws ), where At em is the emission time at the jet 
source frame (e.g., Daigne & Mochkovitch 1998). Also, our 
models are not meant to explain the shape variation of the 
bursts with spectral band (e.g., Norris et al. 1996), which 
would depend on fraction of enegy radiated in gamma rays, 
e. In fact, one expects that e will depend on energy and time, 
as the WS decelerates and the energy is radiated at lower 
spectral bands. 



Finally, the model uses the simplest velocity variation 
which allows an analytic solution: it assumes an instanta- 
neous jump in the velocity of the injected material (step 
function), thus, it makes the simplification that the luminos- 
ity instantly achieves the maximum value. A more realistic 
situation would be a gradual increase in the velocity of the 
injected material which would produce a gradual growth in 
the luminosity. Here we show that the decay of two GRB 
light curves can be fitted by the emission of an decelerated 
or accelerated WS given by these very simple models. Al- 
though, this choice is intended to illustrate the formalism 
it is also true that the analytic solutions allow an explo- 
ration of parameter space and give us an understanding of 
the dependance of the GRB emission on important physical 
parameters like the gamma ratio r and the mass-loss rate 
ratio b. 



6 FRACTION OF THE INJECTED MASS 
LOST BY RADIATION IN THE WS 

In this section, we evaluate the fraction of the mass lost by 
radiation with respect to the mass injected in the WS, Am. 
From equation © the total mass lost by radiation is 



L r {i)c dt ■ 



m(r)7(r)dr - m ws (t)'y ws (t). 



(45) 



In the RHS of this equation, the first term is the total mass 
ejected by the flow that has been incorporated into the WS 
and the second term is the actual mass of the WS. These 
two terms increase with time but their difference remains 
finite because the shocks will weaken with time. 

One can define the fraction of mass lost by radiation as 



Am 



2 dt 



^<ws7ws 



III rh ( T )l{ r ) d ' r 



III m(r)7(r)dr 



(46) 



This fraction Am — > when t — > oo because the total mass 
lost by radiation in the LHS of equation (|45[) is finite. Thus, 
one has to evaluate Am at a finite time to determine the im- 
portance of radiation losses in the dynamics of the working 
surface. 

The total momentum lost by radiation J L r (t)c~ 2 /3 ws dt 
can be obtained from equation ([6]). In particular, in the UR 
regime where /3 WS ~ 1, the mass and momentum losses are 
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0.6 




Figure 6. Fraction of mass lost by radiation as a function of the 
gamma ratio r = 72/71- Each curve is labeled by the different 
values of the mass-loss rate ratio b = 1112/1111. 



the same. In this regime, Am measures the importance of 
both momentum and mass losses. 

If Am <C 1, radiation losses will not change the rela- 
tivistic mass of the WS significantly. Note that in the con- 
stant velocity phase, radiation losses, which change the rel- 
ativistic mass of the WS according to equation ((5|, do not 
affect the velocity /3 WS because the LHS of eq. (|7| is zero 
(since d/? ws /df = 0). Nevertheless, these losses do change 
the momentum of the WS at the critical time and the dy- 
namics of the WS in the second decelerating/accelerating 
phase of Case I and II, respectively. 

We choose to evaluate Am at the critical times 
(eqs. [14] and [26]), which correspond to the end of the 
constant velocity phase. For Case I, T2 C = At, thus, 
J T m(r)7(r)dr = 771171 AraA(l + A); for Case II, n c = 
— Ar, thus, m(r)7(r)dr = mi7iAr(l + A). Then, us- 

ing equations (If 6[) . (|28|) . one can show that in both cases I 
and II, the mass fraction is 



Am 



7„ 3 (6 + aX) 
71 oA(l + A)' 



In the UR regime, where a ~ I , A ~ (far) 1 / 2 , and 7™ 
by equation (|B2|) . this expression reduces to 



Am = I 



r (b + A) 



A(A + r-2) 1/2 (I + A) 



vl/2' 



(47) 



is given 



(48) 



and 



which is a function only of the gamma ratio, r = 72 /71 , 
the mass-loss rate ratio, b = 1712/ rh\. 

Figure [S] shows Am as a function of r for the b values of 
the models presented in Table I. One can see that the mass 
fraction is Am ~ 7% for the model of GRB 070318, but it 
reaches values Am > 50% for the model of GRB 0804I3B. 
In the last case, radiation losses clearly affect the dynamics 
of the WS: mass and momentum are not conserved. 



7 CONCLUSIONS 

We have developed a new formalism that describes the dy- 
namics of an internal WS in a relativistic jet produced by 
variations in the source injection velocity. The WS is formed 
when a fast flow overtakes a previous slower flow. This for- 
malism takes into account that the momentum is not con- 
served because relativistic mass is lost by radiation, in con- 
trast with non relativistic flows. 

Assuming step function variations of the injection ve- 
locity and mass-loss rate we find analytic solutions for the 
WS velocity and luminosity We consider two cases: when 
a pulse of fast material reaches the slow downstream wind 
(Case I); and when a pulse of slow material is pushed by 
fast upstream wind (Case II). In the initial phase, the WS 
is bounded by 2 shocks: one shock incorporates the mate- 
rial from the fast (slow) pulse, the other shock incorporates 
material from the slow (fast) wind. In this phase, the veloc- 
ity of the WS is constant. When the material from the fast 
(slow) pulse is completely incorporated into the WS, only 
one shock remains: in Case I, the WS is decelerated as more 
mass is added to the shock by the slow downstream wind; 
in Case II, the WS accelerates, pushed by the fast upstream 
wind. The WS luminosity in the constant velocity phase is 
constant, and decreases with time when one of the shocks 
disappears. To apply these models to observed GRBs we as- 
sume that a constant fraction e of this energy is emitted in 
gamma rays. 

In the UR limit, the ratio of the Lorentz factors r — 
72/71 and the mass-loss rates b = ma/ 'mi of the relativis- 
tic flows that collide can be obtained directly by fitting 
the light curves of GRBs. As an example, we fit the light 
curves of the GRBs 0804I3B and 070318 with the Case 
I and Case II models, respectively. For GRB 0804I3B we 
obtain the ratios 72/71 = 22.5 and 7712/mi = 233.7; for 
GRB 070318 the fit gives lower ratios, 72/71 = 2.6 and 
mi/thi = 0.1. Since the WS is moving towards the ob- 
server at relativistic speeds (7 ws o ~ 100 — 800), one has to 
correct the observed gamma ray fluxes for Doppler boost- 
ing. Assuming 71 = 100, we estimate mass-loss rates of 
the jets between mi ~ 10 -5 — IMoyr -1 . Note that the 
jet kinetic power is Pkin = ■ymc 2 ~ 10 44_4S ergs _1 . This 
is much smaller than the associated isotropic luminosity, 
Liso = 47rD 2 Fgrb I 'e ~ 10 53_54 ergs _1 , uncorrected for rela- 
tivistic effects. In fact, the isotropic luminosity has no phys- 
ical meaning for our relativistic jet models, which are able 
to produce the Doppler boosted gamma ray flux observed 
at the Earth. 

We also evaluate the fraction of the injected mass lost 
by radiation. For the model of the source GRB 070318 this 
fraction is ~ 7%, while for the source GRB 080413B, one 
finds that more than 50% of the injected mass is lost by 
radiation. Therefore, in the latter source, radiation losses 
change significantly the relativistic mass of the WS and af- 
fect its dynamics. 

The step function variability of the source injection ve- 
locity and mass-loss rate is a simple approximation to the 
real time variability of the injection parameters. Other func- 
tional time variations of these parameters can be easily im- 
plemented in our formalism by integrating the equations in 
§3 numerically. Nevertheless, in the UR regime our analytic 
model is very useful to determine the ratios of important 
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physical parameters (the gamma ratio and the mass-loss rate 
ratio) without introducing any further assumptions. 

In a future work we will study the high energy emission 
produced by these internal shocks, in particular, the fraction 
of energy emitted as gamma rays. 
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APPENDIX A: COLLISION OF RELATIVISTIC 
PARTICLES 

We discuss a simple example of an inelastic collision of 3 rel- 
ativistic particles that shows that when energy is radiated 



in a shock, momentum is not conserved because the rela- 
tivistic mass, that includes the Lorentz factor 7, of internal 
motions, is not conserved. 

Consider two particles in the laboratory with masses 
rrik, velocity Pk, momentum Pk = rrikjkcPk, and energy, 
Ek = mfc7feC 2 , for particles k = 1,2. We consider the col- 
lision between these two particles, where P2 > Pi- In the 
laboratory reference frame the momentum conservation (di- 
vided by c) gives, 



mi7i/3i + m^iPi = (mi + m 2 )7;i27i2/3i 



(Al) 



where 7^2 is the Lorentz factor that corresponds to internal 
motions, and 712 corresponds to the bulk motion of the new 
particle 1-2 with velocity /J12. Energy conservation (divided 
by c 2 ) gives 



mi7i + m 2 72 = (mi + m 2 )7ii27i2- 



(A2) 



The ratio of these two equations gives the velocity of particle 
1-2, 



P12 = 



mi7i/3i + 771272/32 



(A3) 



771171 + 771272 

and the Lorentz factor for internal motions is obtained from 
equation ()A2|) 

771171 + 771272 



7*12 



(mi + 7712)71 



(A4) 



Consider a momentum center reference system where 
the total momentum is zero. This system moves with re- 
spect to the laboratory frame with the speed P12 defined by 
equation (|A3[l . The Lorentz transformation of the energy of 
the two particles before collision in this primed system gives 

E' k = 712 (m fc 7 fc c 2 - c/3i 2 (7n fe 7 fe c/3 fc )) 

= 7n fc 7 fc 7i2 (1 - jflfc/Sia) c 2 , 
and the momentum transformation gives 
Pk = 7n fe 7fe7i2C (p k - P12) ■ 



(A5) 



(A6) 

Using equation (]A3|) . one can check that P[ + P 2 — 0. 

In this system, after the collision, the new particle 1-2 
will be at rest. Before any radiation is emitted, the initial 
energy before collision E' initial is conserved; thus, after col- 
lision it will become the internal energy of the particle 1-2, 
pi 

internal s 



pi 



E'i + E' 2 



(771171 + 777272) c 2 



712 



(mi + m 2 )7ii 2 c 2 = E' in 



(A7) 

where we used equation (|A3|) to calculate E' , and y' il2 is 
associated with the internal motions of particle 1-2 in the 
momentum center reference system. Solving this equation 
for 7- 12 , one gets 



, 77li7i + 777272 

7*ia = 7 1 \ = 7*12, 

(mi + 7772)712 



(A8) 



where we used equation (|A4|) ; i.e., the Lorentz factor associ- 
ated with internal motions is the same in both reference sys- 
tems. Thus, after collision the new particle 1-2 has a larger 
mass, (mi + 7712)7*12- 

Now, let us assume that the particle 1-2 radiates 
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isotropically (in the momentum center system) all the avail- 
able energy, thus, the momentum of particle 1-2 will not 
change, P[ 2 = 0. The maximum energy that can be radi- 
ated by particle 1-2 is 



F' — F' 



712 



[mi7i + m 2 j2 — (mi + 7712)712] , 



(A9) 



where E' = (mi + m^c 2 is the rest energy of particle 1-2, 
and we used equation (|A7[) . 

In the laboratory reference frame, this particle 1-2 that 
has radiated all the available energy now has an energy and 
momentum given by the Lorentz transformations 



E = 712 [E' + c/3i 2 P[ 2 ) = (mi + 7712)7120 



(A10) 
(All) 



Pa = 712 (Pl2 + ~ E o) = ( m l + Tl2)7l2C/3l2- 

Comparing these equations with the RHS of equations (|A1|) 
and (|A2jl one can see that after radiation, particle 1-2 has 
lost energy and momentum. 

On the other hand, the velocity of particle 1-2 after 
radiation in the laboratory reference frame is given by 

= fta. (A12) 



Comparing with equation (|A3[) one can see that, after ra- 
diating all the available energy, particle 1-2 conserves its 
velocity. 

Now, let us consider the collision of three particles with 
rest masses and velocities /3k, respectively, for k = 1, 2, 3. 
We want to obtain the velocity of the particle that results 
from the collision of the three particles. 

As discussed above, /?2 > /3i, such that particle 2 will 
collide with particle 1. The resulting velocity of the com- 
bined particle 1-2 and internal Lorentz factor are given by 
equations (|A3|I and (|A4p . Now we consider the collision of 
particle 1-2 with particle 3 that has 03 > /3i2. 

If the new particle 1-2 collides with particle 3 before it 
radiates away energy, the conservation of momentum gives 

(mi + m 2 )7,i272/3i2 + 771373/33 = 

= (mi + m 2 + m 3 )7ii237i23^i23 , (A13) 
and the conservation of energy gives 

(mi + m 2 )7ii272 + "1373 = (mi + m 2 + 7713)7^123- (A14) 

In this case, where no energy is radiated before the collision 
of the 3 particles, the velocity of the new particle 1-2-3 is 
given by the ratio of equations (|A13|) and (|A14|) 

nad _ (mi + m 2 )7ii272/9i2 + m 3 7 3 /3 3 
P123 — / ; \ ; 

(mi + 7712 )7il272 + 7H373 

_ mi7i/3i + 777272^2 + 771373/33 
TO171 + 771272 + 771373 

where we used equation (|A1[) to obtain the last equality. 
Note that the RHS of the above equation is also obtained 
from the conservation of momentum and energy of the col- 
lision of the 3 particles. If the new particle 1-2-3 radiates 
all its internal energy, as shown in equation (IA12|) . it will 
conserve this velocity. 

We now consider the case when particle 1-2 radiates 



(A15) 




Figure Al. Luminosity jumps J 1 and J 11 as function of the 
mass-loss rate ratio b for two 7 ratios: r — > 1 (solid lines) and 
r = 3 (dashed lines). 

all the available energy before colliding with particle 3. In 
this case, particle 1-2 has the momentum and the energy 
given by equations (|A11 f) and ()A10|) . Then, the equation of 
momentum conservation is 

(mi+m 2 )72/3i2+m373/? 3 = (mi+m 2 +m 3 )7;i237i23/3i23, (A16) 
and the conservation of energy is 

(mi + 7712)72 + 771373 = (mi + m 2 + m 3 )7;i237i23- (A17) 

In this case, when energy is radiated by particle 1-2 before 
collision with particle 3, the velocity of the new particle 1-2-3 
is given by the ratio of equations (IA16|l and (|A17|) . 

a rad (mi + m 2 )72/?i2 + m 3 7 3 /3 3 / A1 o^ 

P123 — -, ; ; ; ■ l A1 °J 

(mi + 7712)72 + 711373 

This velocity is different from the velocity of the adiabatic 
case, equation (|A15[) . Therefore, the collision of particles 
that radiate their internal energy does not conserve mo- 
mentum! This happens because when energy is radiated be- 
fore collision by particle 1-2 above, its momentum decreases 
because its relativistic mass (mi + 1712)7^12 decreases, i.e., 
7»i2 -> 1- 

Therefore, when energy is radiated in a shock in the 
collision of 3 relativistic particles, one cannot obtain the 
velocity of the shock by assuming momentum conservation. 



APPENDIX B: 
CASE 



THE ULTRA-RELATIVISTIC 



Here we will consider the ultra-relativistic (UR) limit, /3 ~ 1, 
where the Lorentz factor 7 > 1. In this limit, a ~ 1, and 
one can expand /3 ~ 1 — l/(27 2 ). Thus, /3 WS — /3i ~ (7ws — 
7i 2 )/27?7ws, and a - 1 ~ ( 7a 2 - 7i 2 )/(27i7 2 2 )- 

We define the normalized Lorentz factors 77 = 7 ws /7i 
and r 2 = r y WB /j2- 
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In the constant velocity phase, the velocity of the WS, 
given by equation (|13|) . can be written as 



/SwsO 



7a + >7? 



(Bl) 



27?7 2 2 (1 + A)' 

where A ~ (for) 1 / 2 (eq. [12] with a = 1). Thus, the constant 
Lorentz factor is 

(l + A) 1/2 7i 7 2 



7ws0 



(B2) 



(7 2 + A 7 ?) 1/2 ' 
and the normalized Lorentz factors are constant and are 
given by 

1/2, 



o_ (l + A)-<-,- 
"i 



and r 2 



(B3) 



(r 2 + A) 1/2 

In the decelerating Case I and the accelerating Case 
II, the normalized Lorentz factors depend on time and are 
simply given by (eqs. [19] and [31] ) 



i 
ri 



T 1 + 2 
T 1 - 2' 



and r| 7 



T 1 



TII+i , (B4) 

where the time functions T l {t) are given by equations (|23p 
and (|35|l . respectively. 

Thus, the luminosity in the constant velocity phase (eq. 
BO]) is 



Lj 



7ujs0 



■7? 



2 7i 2 7Lo 
1 (r? + l) 



(72 — jwso) (7^0 - 71) 



72 
-1) 



71 



r 2 °) 



(B5) 



" 2 7l 2 W) 2 (l + r°) 

where the last term was obtained by solving for A in equation 
(|B3[) . The luminosity of the decelerating Case I (eq. [41]) in 
the UR case simplifies to 

! 3 



L\jR{t) 



(rl + 1) (rl-l) c 



32 



T 1 



(B6) 



47? W) 2 7i ([T 1 ] 2 - 4) 2 ' 

and the luminosity of the accelerating Case II (eq. [43]) sim- 
plifies to 

\3 



A 2 H / + l)(l-r 2 ")^ 



4 72 2 
32A 2 



(r 2 ") 2 



7 2 2 ([T"] 2 -4) 2 ' 
Finally, the dimensional luminosities are 
32T l (t) 



Lttr — I 



([TW-4) 2 ' 

ec 2 rhi/ji for Case I, and I 1 



(B7) 
(B8) 

ec 2 7712/72 for 



where / 7 
Case II. 

We now examine the "jump" in the luminosity at the 
critical time t l c . At this time, ri,j and r 2 ,/i coincide with 
ri,o and r2,o, respectively, given by eqs. (|B3|) . because the 
velocity of the WS is a continuous function of time. Then, 
we define the jump of the luminosities at tf c as 

L'uRitl) _ (r? + l)(l 



J 



1 



r 2 °) 



for Case I, and 



J 11 = 1 



2(r?-r 2 °) 
(r?-l) (1 



r° 2 ) 



2(r?-r 2 °) 



(B9) 



(BIO) 




Figure Bl. Isocontours of the luminosity jump as function of 
the mass-loss rate ratio, b, and the gamma ratio, r. Each contour 
is labeled by the value of the jump for Case I and, in square 
brackets, the value of the jump for Case II. 



Table Bl. Coefficients 



Case 


do 


ai 


02 


bi 


b 2 


I 

11 


-2.0000 
-0.6667 


-1.6724 
-0.5494 


-0.3646 
-0.2651 


0.7635 
0.6247 


0.2431 
0.1768 



for Case II, where we used equations (|B5f) . (|B6|) . and (|B7|) . 
The luminosity jumps are shown in Figure [ATI as a function 
of the mass-loss rate ratio 6 for two values of the gamma 
ratio r: r — > 1 (solid line) and r = 3 (dashed line). One 
can see that small luminosity jumps are predicted for large 
values of b in Case I; the opposite is true for Case II. The 
figure also shows that the summ J 1 + J 11 — 1 because at 
the critical time, the WS of Case I and Case II are each one 
bounded by only one of the shocks that provide the total 
luminosity, L° , in the constant velocity phase. This fact is 
also presented in Figure [BT1 by isocontours of the luminosity 
jumps J 1 and J 11 as functions of the mass-loss rate ratio b 
and gamma ratio r. 



APPENDIX C: THE CRITICAL TIME 
FUNCTION, Tc 

This Appendix describes a procedure to obtain the critical 
time function, 7~ c , from a fit of the GRB light curve. Let 
us assume that the constant velocity phase starts at an ob- 
served time to, then, the decaying wing of the GRB light 
curve starts at t e — to + t c , where t c is the critical time (t e 

7 It can be shown from equations lilt . 1 140 X , H41I I. and (143 b . that 
L° = L 1 + L 11 , is a general result, not only valid for the UR case. 
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measures the end of the constant velocity phase). We further 
assume that the wing of the GRB light curve can be fit by 
a power-law, L — B(t — to)", t > t e . Consider three times: 
t n — nt c + to, n = 1, 2, 3, and the corresponding dimensional 
luminosities, L(T n (t n )), which are given by equation ()B8[) . 
Then, one can construct the relation 



% (r c 2 



Tc (72 - 4) 2 
which can be solved for the time function Tn as 

'T c 2 -4) ? 



-sx, 2 - 



n a T c 



-Tn + 16 = 0, 



2,3. 



(CI) 



(C2) 



These are two quartic equations for time functions T n (T c ,a). 

Equations (|23|) and (|35|) give the relation between the 
time t and T for Cases I and II, respectively. Let us define 
//(T) = T 2 /2 + 4T + 41nT, for Case I, and n"(T) = 
T 2 /2-4T + 41nT, for Case II. Then, one can construct the 
relation 



% = n (T) 



2fj, I ' II (T2) + n I ' I1 (%) = 0. 



(C3) 



When one substitutes T n (T,a) from equation (|C2I) in this 
equation, one obtains a relation H x ' xx (Tc , a) = 0. This re- 
lation gives two distinct curves for the critical time function 
Tc(a) for Case I and II. These curves can be fit by the Pade 
polinomials 



ao + ai (Tc - 2) + a 2 (Tc - 2) z 
1 + 61 (T - 2) + b 2 (T - 2) 2 



(C4) 



where the coefficients are given in Table 2. These curves are 
shown in Figure ICll In particular, the values of the exponent 
a are restricted to the interval [—2.02662, —3/2] for Case I, 
and [—3/2, —2/3] for Case II. Thus, from the slope a, one 
can determine to which case corresponds the observed wing 
of the GRB light curve: Case I (decelerating WS) vs Case 
II (accelerating WS). 



APPENDIX D: DETERMINATION OF 
FUNDAMENTAL PARAMETERS: GAMMA 
RATIO AND MASS-LOSS RATE RATIO IN THE 
ULTRA RELATIVISTIC CASE. 

The gamma ratio and mass-loss rate ratio can be obtained 
by the following procedure. 

I) Given the GRB light curve, one chooses the time for 
the beginning of the velocity pulse, to (this corresponds to 
t = in the models), and the time for the end of the constant 
velocity phase, t e = to + t c , where the critical time is te- 
ll) From a numerical fit of the GRB light curve wing, 
L — B(t — t c ) a , the critical time T for Case I or Case II can 
be obtained from equation (|C3|) . 

Ill) From the observations, the jump J of the luminosi- 
ties at t c , defined by equation (|B9|) and equation (|B10|) for 
Case I and Case II, respectively, are estimated. Also, from 
these equations, writing the gamma ratio r = rj/rij, with 
r? = (T + 2)/(7; - 2) for Case I and r° 2 = (T - 2)/(T c + 2) 



for Case II (see eq. [B4] ). it can be shown that, 

(T + 2) [J (T - 2) - T] 
r (T - 2) [J (T + 2) - Tc] ' 



(Dl) 




-0.5 



Figure CI. Relation between the exponent a and the critical 
time, Tc, for Case I and Case II, obtained from equation dC4ll , 



where this equation applies for both cases I and II. 

IV) Finally, the parameter A = 6 1,/2 r 1/ ' 2 is calculated 
from equation (|B3[) . that is, 



A,= \Tc-J(T-2)Y 
J(l- J)(T c -2) 2 

, _ J (1 - J) (T + 2) 2 



[Tc-J(T + 2)Y 



(D2) 



(D3) 



Therefore, only the ratios A and r can be obtained in- 
dependently by fitting the observations in the UR case. Fi- 
nally, from A and r one can obtain the mass-loss rate ratio 
b = X 2 /r. 



APPENDIX E: DOPPLER BOOSTING 

As found by Lind and Blandford (1985), for a WS in a rela- 
tivistic jet approaching the observer with a velocity /3 WS , the 
observed flux at a given frequency, F v is increased with re- 
spect to an emitter at rest by a factor 8% s , where the Doppler 
factor is <5 wa = 1/[7 WS (1 — /3 ws cos#)], where 8 is the angle 
between the observer and the jet axis. 

In particular, for an optically thin WS, the observed 
flux is given by 



F v 



5 J 

u ws 

LP 



jl'dV, 



(El) 



where D is the distance to the jet source and j' v i is the vol- 
ume emissivity (erg cm _3 s _1 str _1 Hz _1 ) at frequency v' — 
v/5 ws measured in the rest frame of the emitting source, and 
the integration is carried over the volume V' . 

Integrating the flux over a 4-7T solid angle and over the 
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frequency v', one gets 

4n J F v dv' = °^-k j j j' v ,dV'dv'. 

Given the transformation v' = v/Sws, the integral in the 
LHS of the above equation is 

F v dv' = — — / F„dv, (E2) 
o " ws Jo 

therefore, the bolometric flux is 

F„dv = -^L', (E3) 

where L' (ergs -1 ) is the bolometric luminosity or the total 
emitted power of the source. 

In our model, the luminosity of the WS in equation 
(|38[) . L r , is measured in the rest frame of the source of the 
jet. This total emitted power is a Lorentz invariant for any 
source that emits isotropically in its instantaneous frame of 
reference (Rybicki & Lightman 1979). Thus, L r = L' . Fur- 
thermore, since we consider that a fraction e of the bolomet- 
ric luminosity is emitted in gamma rays, Lgrb = tL r (see 
§4), then, the observed flux in gamma rays, Fgrb = e-Fbol, 
is given by 

Fgrb = ^ W (E4) 

This equation can also be written in terms of the non di- 
mensional luminosity in equation (|39|l as 

Fgrb = L e rh^c . (E5) 

The case of an optically thick WS, discussed by Lind 
& Blandford (1985), follows in a straight forward fashion. 
One also obtains the 5^, a factor between the emitted and 
observed bolometric fluxes because the extra <5 WS factor is 
due to the integration in frequency. 
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